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ABSTRACT 


Local and average heat transfer coefficients were deter-~ 
Meem@ed £LOr a right circular cylinder in an oscillating flow. 
Spanwise platinum heater strips were used to heat the cylinder 
isothermally over the lower 180 degrees from the front to rear 
stagnation point. The four inch diameter cylinder was posi-~ 
tioned both normal to and at 45 degrees to the flow direction. 

Data was gathered for diameter Reynolds numbers from 
100,000 to 300,000. Large amplitude oscillations were imposed 
upon the mean flow using a rotating shutter arrangement. Fre- 
mremeres Of 0, 5, 10, 22, 50, 100 and 126 Hertz were investigated. 

wMorenormal Elows, iocal heat transfer coefficients in the 
wake and average Nusselt numbers were enhanced above values 
Memecsteady flow, for runs having large amplitude. In the 45° 


flows, no significant change in heat transfer was noted with 


Saeiillating flow. 
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Pee OUUC TION 


A. NORMAL CYLINDER 

Mereed comvective heat transfer from a cylinder in uniform 
flow has been a topic of interest to researchers since the turn 
of the century. In 1914 Professor L.V. King presented his 
relationship to describe the heat transfer from a cylinder in 


laminar crossflow as a function of air velocity: 


MeEnoOwetweslighntiy modified by Collis and Williams in 1959, 
it still serves today and is the basic equation of hot wire 
anemometery. Numerous other correlations are available and 
Mmemasstandard part of the heat transfer literature [Ref. 1]. 

A topic of increasingly greater interest to researchers 
and designers is the effect of an oscillating flow on forced 
convection. Gas turbine designers are intensely interested 
in upgrading the performance of gas turbine blades through 
materials research and blade cooling. Gas turbine blades 
experience a hot, oscillating flow environment as they pass 
through the wakes of upstream blades. 

Heat exchange designers involved in waste heat recovery 
from diesel engine exhausts are also concerned with the effect 


of flow oscillations on the exchange of heat to and from 


cylinders. 
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Piaaalomdataiacmpeen made available to date for a cylinder 
in oscillating flow. Data is available for heat transfer from 
meblate plate in oscillating flow [Refs. 1, 2 and 3] and fot 
metic toltl in oscillating flow [Ref. 4]. Data for a cylinder 
mother than uniform flow is generally for a vibrating cylinder 
Mme Or without a net £Llow [Ref. 5] or for a cylinder with 
screen induced turbulence of less than about ten percent [Refs. 
me, & and 9]. This study concerns itself with the effect of 
oscillations on the order of fifty percent of the freestream 
velocity on local heat transfer rates. 

Yawed cylinders are of interest to aerodynamicists because 
of the correlation of a yawed cylinder to a swept back aircraft 
wing. Researchers and designers may also encounter yawed cylinders 
in heat exchangers and other applications. 

No information is currently available on the effect of 
[ow Oscillations on a yawed cylinder. Viteri [Ref. 10} deter- 
Mined the average heat transfer coefficients of yawed cylinders 
Mmimiiitorni cLlow, and Kraabel [Ref. 11] determined the local 
heat transfer coefficients around yawed cylinders in uniform 
flow. These papers provide a basis for validation of the yawed 
cylinder results reported here. 

It has been demonstrated by Miller [Ref. 2] that augmented 
convective heat transfer can be attributed to the reduction of 
laminar to turbulent transition Reynolds numbers. Despard 
[Ref. 12] showed that freestream oscillations promote laminar 


boundary layer separation in the presence of an adverse pressure 


al 





eo Web annmine s|Rero £5] produced data that revealed that 
the starting vortex associated with periodic separation from 
an airfoil in oscillating flow produced an unexpectedly strong 
boundary layer reattachment which markedly improved the aero- 
dynamic performance of the foil. 

In summary, the effects of early laminar to turbulent boun- 
dary layer transition and the destabilized boundary layer pro- 
duced slightly augmented heat transfer while periodic vortices 
shed from the separation point can cause periodic reattachment 
of the hydrodynamic boundary layer. This can be expected to 
greatly enhance the heat transfer in the afterbody region. 

The purpose of this investigation is to determine the effect 
of large flow oscillations on the local heat transfer coeffi- 
cients around both normal and yawed cylinders. Additionally it 
is desired to develop a correlation formula for the average 
Nusselt number in oscillating flow that includes the effect 


of yaw angle, Reynolds number, and flow frequency and amplitude. 
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PDO webeoewee ELON OF @ePPARATUS 


A. MODEL 

The model used in the investigation was a four inch dia- 
meter cylinder. It was fabricated from one half inch wall 
thickness fiberglass-phenolic tubing two feet long. The 12 
inch mid-span section was instrumented with 16 platinum heater/ 
resistance thermometers oriented parallel to the axis of the 
cylinder and spaced 12 degrees apart from center-to-center. A 
sketch of the model showing the orientation of the heater strips 
meee Shown in Figure 1l. 

The platinum strips used were 12 inches long by one-quarter 
Of an inch wide by one-thousandth of an inch thick. Each strip 
subtended an arc of 7.16 degrees of the circumference of the 
cylinder. The strips were spot welded to brass studs that were 
machined flush with the cylinder on the outside and protruded 
into the interior of the model for attachment of the voltage 
sensing and current carrying wires. Figure 2 shows a cross 
section of the cylinder wall through the studs and strips. 
Voltage drop across the platinum strip was measured using the 
voltage sensing wires attached to each stud. Platinum was 
used to both heat the cylinder and determine the surface tem- 
perature of the cylinder because of its large linear resistivity- 
temperature relationship. The resistance of the strips was 


determined from the voltage drop across the strip which was 
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Bem to be used to find the temperature. The amount of heat 
Semvected at “each Location was proportional to the power dis- 
Spaced bpy the Strip which was equal to the product of the cur- 
rent squared and the strip resistance. 

The entire outside of the cylinder, including the platinum 
strips, was covered with a thin layer of hard epoxy. The hard 
epoxy layer was machined and polished on a lathe to reduce the 
thickness over the strips to less than one-hundredth of an inch 
and to produce a smooth heat transfer surface. A photograph 
of the completed model showing the mounting stud and bracket 
and electrical leads is shown in Figure 3. 

Two copper-constantan thermocouples were mounted on the 
inside surface of the cylinder at 90 degrees and 2/70 degrees 
from the front stagnation point. These thermocouples were used 
to monitor the interior wall temperature of the cylinder in order 


to estimate the heat loss to the inside of the cylinder. 


B. WIND TUNNEL 

The experiment was carried out in the oscillating flow wind 
tunnel of the Naval Poe dae School. Photographs of the 
Facility are shown in Figures 4 and 5. The oscillating flow 
wind tunnel is a low speed, open circuit wind tunnel with the 
inlet bell extending out of the building. The tunnel inlet 
bell is eight feet square and the test section is two feet 
SteucrmeGoOvratne a lo.) Contraction ratio. Three high solidity 


Seteens Located in the inlet section upstream of the nozzle 
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megduee ftreestream turbulence intensities of less than 0.5 
percent of the freestream velocity for uniform flow. 

The wind tunnel drive consists of two axial flow fans in 
Series, each of which has an internal, 100 horsepower, direct 
connected, 1750 RPM motor. The fan blades are internally adjust- 
able through a pitch range of 25 to 55 degrees, providing a 
wide operating base. A set of variable inlet vanes, located 
immediately upstream of each fan, are externally adjusted to 
Meovide control of test section velocity. These vanes, of radial 
@emeiguration, preswirl the air in the direction of fan rotation 
to reduce fan capacity. The total range of frseestream velo- 
cities possible is from ten to 250 feet per second. 

A plan view sketch of the wind tunnel showing the location 
of the fans and inlet vanes in relation to the test section is 
shown in Figure 6. 

Two fundamental methods are available to superimpose oscil- 
tating flow, on a mean velocity. Both Nickerson [14] and 
Hori [15] introduced oscillations by oscillating their models 
in a steady flow environment. This method severely restricts 
the range of attainable frequencies because of mechanical com- 
plications and instrumentation difficulties. 

The other approach is to provide an oscillating flow environ- 
ment over a stationary model. Hill [16] used a sliding shutter 
to impose oscillations on the freestream but was limited to 
low frequencies because of mechanical difficulties. Feiler 


Boa, Yeager [17] used a siren mounted upstream of an eight inch 
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diameter wind tunnel to expose small models to frequencies 
ranging from 34 to 300 Hz. with root-mean-squared amplitudes 

of up to 65 percent of the mean velocity. The most successful 
method of obtaining oscillating flows with a large range of 
Frequencies, amplitudes and flow velocities in a moderate sized 
tunnel is that employed by Karlson [18] and later by Miller 
[19]. A rotating shutter valve, immediately downstream of the 
test section is used to superimpose a periodic variation of 
velocity on the mean flow. The method used in this investiga- 
meonm 1s identical to that used by Miller. 

The shutters consist of four horizontal, equally spaced, 
shafts mounted in a plane just downstream of the test section. 
The shafts are slotted to receive blades of various widths, 
forming a set of four butterfly valves spanning the wind tunnel. 
Figure 7 is a picture of the shutter valve assembly. The 
Shutters are driven synchronously by a pulley and belt arrange- 
ment and a vari-drive motor and gearcase, and can be varied 
meome0.1 — 250 Hz. The amplitude of the oscillation is a func- 
tion of the width of shutter blades, frequency of rotation of 
the shutters and the freestream velocity. Five sets of blades 
are available producing tunnel blockages of 33, 50, 66, 83 
and 100 percent of the flow area in the fully closed position. 
In this investigation blades providing 66 and 100 percent 
blockage were used resulting in amplitudes of 5 to 69 percent 


of the mean freestream velocity. 
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GiemibesNSectLiOn Of Ehe tunnel is fabricated from two-two 
mienethitck, 25 inch wide by 18 feet long pieces of aluminum 
as the upper and lower walls and three-two inch thick stress 
relieved lucite panels on each side as side walls. The front 
panels are hinged at the top and can be hydraulically opened 
and closed for access to the test section. 

For this investigation the rear center, lucite panel has 
been replaced by a two inch thick plywood panel to facilitate 
mounting the model. The model was mounted by a chordwise 
extension held in a rotatable clamp mounted flush in the center 
of the plywood panel. The other end of the model is supported 
Meee o/s inch thick “Y" shaped, faired bracket, that is affixed 
to the tunnel ceiling and floor, and rests against the front 
ieee tte panel. The model is shown in place with the “Y" bracket 
clearly visible in Figure 8. Electrical leads were taken out 
through the plywood back panel where they were connected to the 
power supplies and the voltage drop sensing circuits. The 
model was yawed by affixing extension pieces to each end that 
were angled and sized to completely span the tunnel at each 
angle of yaw. 

Power to heat the strips was controlled by sixteen variacs 
that powered 16 individual DC power supplies (see working cir- 
GCuit, Figure 9). <A precision shunt (R=0.001 ) was placed in 
Series with each strip and the voltage drop across the known 
resistance was used to determine the current. The voltage 


drop across each strip was measured using the potential sensing 
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leads that were attached to the studs on the inside of the model. 
Instantaneous velocity measurements for the purpose of free- 
stream turbulence determination in uniform flows and amplitude 
and frequency measurement in oscillating flows were made using 

a linearized constant temperature hot-wire anemometer described 
in Reference 20. An oscilloscope and a Ballantine True R.M.S. 
meter were used to monitor hot-wire signals. Mean velocity of 
the free stream was measured using a Prandtl Probe and a water- 
filled micromanometer. 

Figure 10 shows a pair a typical velocity wave forms. The 
vertical axis is velocity and the horizontal axis is time. The 
roughness of the upper tracing is typical of low amplitude oscil- 
lations which are significantly more difficult to produce in 
true sinusodial form. 

The freestream temperature and the model internal wall 
temperature were measured using Beane aoa s atin thermocouples 
and a Newport digital temperature indicator. 

The frequency of oscillation was also monitored with an 
optical chopper driven by a shutter valve shaft and a digital 


Geuneter. 
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ie EXPERIMENTAL PROCEDURE 


fee CALIBRATION OF MODEL 

The platinum strips were calibrated as resistance thermo- 
meters by determining their electrical resistance as a function 
of temperature in an electric oven. Seven separate tempera- 
tures were used, evenly spaced between ambient (60°F) and 
200°F. Although the platinum heater-thermometer strips were 
extremely pure and are expected to have a linear resistance- 
temperature characteristic, calibration was carried out to 
determine the exact slope and intercept of each strip. 

The calibration circuit employed (see Figure 9) was powered 
by a low impedance DC voltmeter calibration standard. This 
proved to be an extremely precise, drift free power supply. The 
voltage drop across the precision resistance (1.0 gi OE A mwas 
recorded as was the voltage drop across the heater strip itself. 
The current was then computed from the voltage drop across the 
precision resistance divided by one ohm, and the resistance of 
the strip for that temperature was computed from the voltage 
drop across the strip divided by the current. Temperatures 
in the oven were monitored with a Newport Laboratories Digital 
Temperature Indicator and a copper-constantan thermocouple. 

After determining the resistances as a function of the 
temperature a straight line was fitted to the data using linear 


regression. The resulting equations were used throughout the 
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PeseduetwLecxperimentation and data reduction to determine sur- 


face temperatures. 


Bee ADJUSTMENT OF SURFACE HEATERS 

Before each series of runs the heater strips were energized 
and the wind tunnel was run with a low velocity uniform flow 
in order to allow the model substrate to come to thermal equili- 
brium. Once this was accomplished the shutters were turned on 
and adjusted to the desired frequency and the velocity was 
adjusted to the magnitude corresponding to the Reynolds number 
desired. At this point the model was again allowed to come to 
thermal equilibrium. 

Once this initial thermal equilibrium was achieved the 
voltage drop and current through each strip was recorded and 
inserted into a predictor-corrector program (see Appendix F), 
This program calculated a new required voltage drop based on 
desired strip temperature, current strip temperature, and the 
known resistance-temperature calibration curve for that strip; 
the strip was then adjusted to the new desired voltage drop 
through adjustment of its power supply. Following adjustment 
of the strips the model was again allowed to return to thermal 
Serb tbpDrium. 

These adjustments were performed as many times as necessary 
until all the strips were at the same temperature, ensuring an 
isothermal heat transfer surface that minimized heat conduction 


between strips. 
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The temperature for a particular run was determined from 
a sampling of the strip temperatures before the first set of 
adjustments was made and was selected to minimize the number 
of adjustments required. In every case however the temperature 
was chosen to be over 110°F in order to maintain at least a 
sor temperature difference between the strip and the ambient 


temperature. 


freee REESTREAM CONDITIONS 

A series of runs was made with both the 45° yawed model 
and the normal model covering a wide range of operating condi- 
tions. Each model was tested with both the six inch blades 
(100 percent tunnel blockage) and the four inch blades (66 per- 
Seme tunnel blockage), resulting in oscillation amplitudes from 
Five to 69 percent of freestream mean velocity. Each combina- 
tion of frequency and amplitude was run at two freestream velo- 
cities corresponding to nominal diametral Remolds numbers of 
bow o00 and 150,000. Frequencies investigated were 0, 5, 10, 


meee, 100 and 125 Hz. 


Dee DATA REDUCTION 

The data reduction technique is described in Appendix A 
and a sample calculation is presented in Appendix B. A detailed 
energy balance is necessary in order to determine what portion 
of the total Joulian heat generated by platinum resistance 
heater is transferred from the surface by convection. Part of 


the energy was conducted through the thin epoxy coating and 
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convected directly to the freestream or was radiated into the 
environment. The balance of the energy was conducted into the 
substrate of the cylinder where some was conducted and con- 
vected into the center of the cylinder, some was lost to the 

ends of the cylinder, and some was conducted to the gaps between 
the strips where it was convected to the freestream or was radia- 
ted to the environment. 

The energy lost from the substrate side of the surface 
heaters is hard to determine analytically since it consists of 
both conduction to the intraheater gaps and subsequent loss by 
convection and radiation as well as conduction to the inner 
surface of the model. Therefore a two-dimensional Teledeltos 
paper model was constructed to determine the relationship 
between the thermal resistances to the inner and outer surfaces. 
This procedure and the results are discussed in Appendix D. 

The results of the Teledeltos paper experiment were placed in 
equation form and solved interatively for the external covec- 
Meo heat transfer coefficient, h, for each strip. 

The other major loss of energy from the system is by radia- 
tion to the environment. These losses can be appreciable and 
an accounting for this loss was developed as part of the data 


reduction program. Details are described in Appendix C. 
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Pee ooo AND DISCUSSION 


A. STEADY FLOW 
Six steady flow runs were made, two with the model normal 


to the flow at Ne = 150,000 and 300,000 and four runs with 


the model at 45° to the flow at [Pu 000n. 190,000, 150,000, 


and 300,000. The results of these runs are shown graphically 
in Figures 13 through 18. Results of the model at 90° (normal 


Bemtene flow) at Npo = !507000 > are compared ins Figure 19 with 


the theory of Froessling [Ref. 24] and the experimental 


Poeewmeits of Schmidt [Ref. 25] for Neo = 170,000. The discrepancy 


between theory and both of the experimental results is expected 
and is a consequence of the effect of freestream turbulence 
near the stagnation point [Ref. 26]. 

Comparison of the results of the present experiment with 
those of Schmidt indicates two interesting points. The two 
agree within 20 percent up to an angle of 150° from the 
stagnation point, beyond which they diverge markedly, differing 
by 70 percent at 180°. This sharp increase in heat transfer 
above that noted by Schmidt is typical of all the results 
observed in the present investigation and shows no correlation 
with frequency number or Reynolds number. Agreement within 
20 percent is good however and this lends a high degree of 


confidence in the experimental techniques and results. 


ae 
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Nomlocal heat transfer data is available for the steady 
flow case with 45° yaw. However, the data of Viteri [Ref. 10] 
can be interpolated to arrive at an average Nusselt number of 
128.3 at a Reynolds number of 31,900. At a Reynolds number of 
35,630 the present experiment produced an average Nusselt number 
of 122.1. This close agreement lends confidence in at least 


the average Nusselt Numbers obtained in the A yawed case. 


Bee NONSTEADY FLOW 

Two thirds of the nonsteady freestream conditions did not 
produce substantial alteration in heat transfer when compared 
to the steady flow case. Figures 20 through 25 depict the 
cases with the 90° model that show little or no altered heat 
transfer, and Figures 26 through 39 depict those of the a 
model. In each of these figures the steady flow data has been 
plotted along with the nonsteady data for comparison. 

Twelve of the 38 freestream conditions investigated pro- 
duced substantially altered local heat transfer rates. These 
results are shown graphically in Figures 40 through 49 for the 
90° model, and Figures 50 and 51 for the 45° model. The 
changes in local Nusselt number took place almost exclusively 
in the afterbody region of the model, and resulted from altered 
separation and flow patterns in the wake of the cylinder. Fig- 
ure 40 is a typical example of these results and demonstrates 
enhanced heat transfer in the wake area of the 90° model. The 


data reported in Figures 48 through 51 exhibit diminished heat 


Pmegeter rates. Figure 50 is typical of this class of data. 
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MmoGdeaetoesim=iplainy the analysis of the results, the 
average Nusselt number was calculated for each freestream con- 
dition (see Table I). Figures 52 and 53 report average Nusselt 
numbers versus Reynolds number for the 90° model and the 45° 
model respectively. In the 90° case CFiewre 52) .a curve has 

O57 6 


been drawn representing Nu = ‘Re and the Amplitude and 


Storuhal numbers have been annotated for each data point. Of 
interest is the fact that data with higher Amplitude numbers 
are those which tend to exhibit augmented average heat transfer 
rates. Moreover, higher Amplitude Numbers seem to correspond 
to greater heat transfer augmentation. Although the data is 
Sparse, a pattern of larger Nusselt Numbers for larger Ampli- 
tude numbers at the same Reynolds number seems to exist. Three 
apparent anomalies exist associated with the data from run num- 
Mmems 15, 1/7 and 20. It is possible that an error in recording 
the amplitude data was made and that the Amplitude numbers for 
runs 17 and 20 were reversed. If this were true then both 

runs 17 and 20 would conform to the apparent pattern. When the 
Frequency numbers, Nusselt numbers and Reynolds numbers of runs 
15 and 17 are compared, it is seen that they are very similar 
Peewthat the Amplitude numbers are considerably different. If 
the previously noted pattern for higher Nusselt number for 
higher Amplitude number is true then the Amplitude number of 


run 15 is suspect and may be too low. 


ey 
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Average Nusselt number data taken for the 45° model is 
shown in Figure 53. No apparent pattern of heat transfer altera- 
tion is evident in this plot even though Amplitude numbers as 


high as 0.44 and 0.65 appear in the data. 


Gee NUSSELT NUMBER CORRELATION 

The data obtained in this investigation is too sparse to 
allow a correlation between heat transfer rate and freestream 
conditions to emerge. Only a possible pattern is suggested by 
the results and more data covering a wide range of Reynolds 
numbers and Amplitude numbers is necessary to produce such a 


@oeerelation. 
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V. CONCLUSIONS 


Pibenouehn there 1S ansurficient data to draw any definite 
conclusions or to derive a correlation, several trends may be 
observed. .In the 90° case, significant increases in the local 
Nusselt numbers are seen in the wake region and are due to 
changes in the separation and wake structure. Moreover, aug- 
mentation of local Nusselt Numbers seems to increase with 
increasing frequency as suggested by the data shown in Figures 
mueethnrough 47. Increases in average Nusselt Numbers is appar- 
ently related to the Amplitude numbers as discussed above and 
Sumewm in Figure 52. With the exception of the three points 
associated with runs 15, 1/7 and 20, discussed in Chapter IV, 
the pattern exhibited is for higher average Nusselt numbers to 
be associated with greater amplitudes. 

In the case of the 45° yawed cylinder, no general patterns 
are discernable and it appears that local Nusselt numbers 
exhibit no significant change over a wide range of frequencies 
and amplitudes. 

This investigation has served to initiate a study of heat 
Beemarer £rom a cylinder in oscillating flow. Additionally 
it has validated the techniques to be used for further study. 
More extensive data collection will be required before defini- 


meavemconciusztonms can be drawn or correlations made. 
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ARE SBE. 


DATA REDUCTION 


The data was reduced to the form shown in Table 2 by a 
specially written Fortran computer program on an IBM 360 conm- 
puter at the Naval Postgraduate School. 

The raw data consisted of pitot-static tube AP (inches 


H,0), freestream temperature, T Cae internal model tempera- 


KF? 
tures (heater side and "other" side, °F), frequency of oscil- 
iemeon (Hz), amplitude ratio (4) and the voltage drop across 


(volts) and current passing through (amperes) each strip. The 


total power dissipated by each strip is simply: 
m5 4a xe x OE 


where 3.4149 is a conversion factor from watts to BTU/Hr. 
From the measured current (I) and voltage drop (E) each 

individual strip temperature is determined by comparing its 

resistance to its known temperature versus resistance charac- 


teristics. The resistance of a strip is: 
E 
R= 
I 


From the linear temperature-resistance relationship obtained 


during calibration, the strip temperature can be obtained: 
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where C and S are the known ordinate intercept and slope 
obtained from the individual resistance versus temperature 
eabibration curve. 

The other unknown temperature required by the data reduc- 
tion program is the inner wall temperature, TT: The inner wall 
temperatures are known at two points. Thermocouples are mounted 
on the wall at locations 90° and 270° from the Feat nat. On  pOlnEe 
at mid-span (see Figure 2). The inner wall temperature is 
assumed to vary linearly between these two measured values. 

With this information h, the heat transfer coefficient may 
be calculated. The following formula (developed in Appendix 


Dpmets used: 


h is contained in the expressions for Ry» RO/Ry> and 


R,/R 


Mearatively to obtain h. 


p 2s 1s shown in Appendix D. This equation was solved 


With the local h for each strip, the local Nusselt number 


momecatculated from: 


_ bD 
Nu = k 
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APPENDIX B 


SAMPLE CALCULATION 


HOmmeEne DUrpOSe Of this sample calculation strip Nr 6 from 


heme run Nr 20 was selected. 


The information recorded is: 





Voltage drop Beam lO 2 VOLES 
Current I = 7.47 Amps 
Pitot-Static Pressure HW = 4.16 inches HO 
Ambient Temperature TA = 64.0°F 
Inside Temperature: 
Heater Side TIHTR = 105.0°F 
Other Side POP = 7 7 OSE 
Frequency Freq = 50.0 CPS 
Amplitude Ratio Ny = 0.31 
Blades Medium (4") 
Model 90° to stream 
Barometric Pressure 30.04 inches Hg 
Freestream Mean Velocity 
2 
_ Jwater, _ 
V | 2g au 1» | 
Ne 
62.4 PS 
= |2(32.2 =£) (4.16 in H,0) (75+44) t- - 1 
S " 0.07654 —— 
Et 
eee 
sec 


SA 





Freestream Reynolds number: 
fete 


, _ vp _ eS ce as Cores 3 rt) 
Re Vv 0.000166 ft2/s 
= 269,800 


Temperature of the strip: 





= - 
dh 
S 
E 1.622 Volts _ 
1 = ee = 0.21/71 ohms 


Cea early G95 5007 


Se Oe OUOSIUGZ5 


Therefore, 


Pee ee Oe liao oO Oo 
i 0.000390625 eee 
Total Heat Dissipated by the strip: 
G. = EL = (3.4149 eee SSN | O72 Volts) (7.47 Amps) 
Ss Volt-Amp-hr ‘ : 
= 41.38 Btu 
Also: 
R Tex. 
eee _ _D _D 
Gaose cr, TR CL + + (Ts - TLR 7 
D O ale 
ko eh 
oS ee 
D kh,A 
R 
0 
Re 5c ieor 


53 





After substitution, the equation may be solved for Den: 
(the constants ays aos by, b, and b4 are given in Appendix D). 


By iteration: 


Dicit 


Zs O 


h. = 36.9 -—-_"— 
hr - ft - “F 


E 


The radiant heat transfer coefficient: 
ee eee Ces 7) 
rE S A S A 


Assuming T, is the same as the fluid temperature: 


A 


oe Gone cms Cs ONS 4 524.0) 


i 
il 


x (570.5° + 524.0°) 


Btu 
hr - ft - F 


0.824 


Therefore, the local convective heat transfer coefficient is: 


S69 = 02624 


Btu 


hr - ee - OF 


Sire 
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tme tocal Nusselt number: 


hd ese. 19100. 333) 


ee ae 0.014738 


G59 


a 





EG ir 


temperature T, 


the net 


where 4 is the Stefan-Boltzmann constant, 


anwea of 


The 


vection 


APPENDIX C 


RADIATION 


PeUdihormlymconvex otrdy body with emissivityG@, at a 
radiating into an enclosure of temperature T,, 


radiant heat exchange is: 


i 


dha son = To 


and A is the surface 
the emitter. 
above equation can be put in the same form as the con- 


equation: 


qd. = h A(T - De) 
where h. [iMac oni Octive Tadi ation Neat transfer coefficient 
paven by: 
ho = @s (7 + T_)(T* + T.*) 
r a 


Daring 6ee=- 0.72 £OE epoxy coating over aluminum [Ref. 22] 


and the following typical temperatures: 


Ties ieO © = 580° R 
T= 60 F = 520°R 
Gives: 
— =3 2 2 
Orme 14) CUGl (580 + 520)(580, + 520°) 
= 0.835 a 
hme £e sk 
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iittomvaluemraArrhs Within two to 1O percent of the total 
heat transfer coefficient and therefore cannot be neglected. 
Therefore, the convective heat transfer coefficient, h, is 
Obtained by subtracting the radiant heat transfer coefficient, 


ho» from the experimentally derived total heat transfer coef- 


Mrerent, that is, 
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APPENDIX D 


CONDUCTION LOSSES 


Although the thermal resistance of the fiberglass substrate 
is very high compared to the thermal resistance between the 
Seep and the fiuid, some heat is lost to the interior of the 
model. This heat loss must be deducted from the total heat 
mmm iat Order to arrive at the actual heat transferred to the 
meee neat Loss due to conduction in the longitudinal axis direc- 
tion and radially through the small brass studs was neglected. 

A Teledeltos paper model was used to estimate the heat loss 
to the interior of the model. Symmetry was used to reduce the 
size of the required model (see Figure 11). 

Teledeltos paper is a graphite covered paper having a nomi- 
nal resistivity of 1000 ohms per inch square and is highly 
isotropic. Teledeltos paper is extremely useful for two- 
dimensional heat conduction modeling in cases such as this 
which have boundary conditions and a physical shape that are 
mreemvetiient for an exact solution. 

The heater strip itself is simulated with a strip of silver 
conductive paint that is placed on the paper in exactly the 
Same position as the platinum strip is on the actual model. 

The inside of the model segment is considered isothermal so 


the analogous surface of the Teledeltos paper model was 
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fofeaeetnitorm potential with a conductive strip of silver 
meee tne Left and right hand adiabatic boundaries, between 
strips, are easily modeled by leaving them bare. The convective 
boundary at the gap between strips is modeled by allowing the 
Teledeltos paper to extend beyond the actual cylinder surface 
for a distance to simulate the convection resistance. The 


required distance is determined by: 


where k is the conductivity of the actual substrate, h, is the 
convective coefficient being modeled and S is the scale factor, 
in this case, S = 48. Additionally the excess length was slit 
into eight equal width strips lengthwise to force the flux 
lines parallel to each other and perpendicular to the surface. 
In order to model h's of 10, 20, 50 and 100 Btu/hr-ft*-F, L's 
of 14.4, 7.2, 2.88 and 1.44 inches respectively were used. 

A strip temperature of 120°F with an inner temperature 
of 100°F and a fluid temperature of 60°F were used. These 
temperatures were scaled to appropriate voltages that would 
avoid excessive Joulian heating of the Teledeltos paper. 

In operation the voltages were measured at 1/7 equally 
Spaced intervals along the line corresponding to the cylinder 
surface and at five equally spaced points along a linecor- 
responding to 0.08 inches from the inner surface. This allows 


Pi m—lciecniiiatlone or the Fatio of heat loss to the inside of 


Sy, 





Ehe model, compared to the heat loss indirectly through the 
miei vyia the substrate). 

Considering the lumped electrical resistance network 
analogous to the heat transfer problem (see Figure 12) it is 
noted that the power input to the strip is equal to the sum of 
the losses to the air, to the inside of the cylinder to the 


gap between the strips, or: 


De ned. 4 
Since in general: 
ee 
Ro - 
this becomes: 
To = Ty T = Ty To = Ty 
= OR awa) "a 
: D 0 I 
Where To Ty and Ty are the temperatures of the strip, the 
fluid and the inside respectively and R ReweandeR. are the 


Dew 0 il 


thermal resistances from the strip to the air, from nena 
Emeeomeh the gap to the air (Ro = Ro + Ri) arcmin Olle | Eire 
Substrate to the inside of the cylinder respectively (see 
Figure 12). 


This equation can be manipulated into the following form: 
eee Rp Boer oD 
Pinel § pei * Rg | Ro 
D 0 D i 
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Mimemresistance Eo heat loss from the strip directly to the 
gar , Ry» Peto OmENe ME rnm bayer of epoxy on the strip and 
the convection film resistance. The conduction resistance 
due to the epoxy is equal to £/kA where £ is the thickness of 
the epoxy layer, k is the thermal conductivity of the epoxy 


and A is the surface area of the strip. The convective resis- 


tance is equal to l/hA, where h is the convection heat transfer 


coefficient. Summing these we obtain: 
5 eet 
D khA 


ime quantities RO/Ry and RL/Ry are derived experimentally 
From the Teledeltos paper analogue and are in the form of 


polynomial equations in h: 


R 
= = 1.546744 + 0.0320672 h 
D 
Re 2 3 
= = 0.106105h + 0.0047559h~ + 0.0000079h 
D 


These resistances were used in the heat transfer equation 


to solve for h (see Appendix B). 
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APPENDIX E 


UNCERTAINTY ANALYSIS 


For the purpose of this analysis, the method of Kline and 
Meelintock [Ref. 23] will be used, assuming that all variables 
have a normal Gaussian distribution. Possible dimensional error 
in strip size, placement and other dimensions were considered 
so small as to be negligible. The tabulated values of k and 
e were assumed to be correct. All temperatures measured with 
thermocouples were assumed to have an uncertainty, due to 
readout error, of 0.1 degree F. 

The strip temperature uncertainties derive from the 
70.003 mV uncertainty in reading the digital voltmeter used in 
@aerbration and the +0.002 V and +0.01 amp uncertainty in 
reading the meters during data runs. | 


The uncertainty in measured resistances during calibration 


1/2 
Z Z 
R a + (=) 


+0.00031 ohms 


is 


> 
bd 
ll 


imemresisctances obtained in calibration are used to obtain 


a straight line equation for temperature vs resistance. 


Rese oe beak <0 
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The uncertainty 


Re The uncertainty 


a0S 


Since the strip 


The uncertainty 


4 
Sots | 


AT. = 


The values of Ro 


model. 


any two points on the Teledeltos paper was +0.02 volts. 


. O 
Secesponas to a uncertainty in temperature of +0.2 . 


Rea R 


o/ Rp 18 given ie 


+0.54°F 


OamGmebomossumead to be of the same order as 


i 


temperature during a data run is given by: 


Of oss 
R - -AT 
s |} (42) + (7) 


+0.00000064 


GOietvis given by: 


eye 


2 2 Ze 
Se aos, a (es 
Be I S 


4- 


1/2 
ue Z 
T + (=) | 


/Ry were obtained from a Teledeltos paper 


Uncertainties in measuring a potential difference between 


as 


For 
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mrere ADL 1s the direct convection surface area times (T, - T.) 
mmemonol is the indirect convection area times (IT local - T._). 
Mi local is determined from the Teledeltos paper model.) 


The uncertainty of SADT is given by: 


SADT = SADT (AAP) 
= +0.0011 
memilkarly for ADT: 
A ADT = ADT (ae 
= +0.0021 
Now 
| ye 
Aa Bo 7 a0 (ADT) i (ASADT)* 
Ry R ADT SADT 
= +0.0088 
RL/Ry heme inven py : 
Ry LADT k (Tt - T.) 


The uncertainty of IADT is given by: 


ALIADT AAD ANS) 


il 


IADT (= 


+0.032 
44 


4 
/ 
{ 

i 





ihnus the uncertainty of RI/Rp imeeeiven by: 





- 1/2 
R R D 2 2 
Pelee AND TE ATcAD T AAT ea 
4 = aa, t Sr + To = err 


Ps) 
Gs) 


fers given by: 


Q 
i 
es] 
KH 


The uncertainty in q is given by: 





ny 2 
2 Zz 
Aq = q A a (45) | 
Sr OnOG 
a ee 


Because of the complexity of the formula used to determine h, 


imgdas tound by iteration of the following implicit equation: 


hk i i 
ea) 
kt+hf S iE ee eet S I Blh + B2h- A B3h> 
The uncertainty of h was determined by: 
2, 2 2 
AT AT AT 
San B2h 2 
Ah = [¢ or) + 6 a7) Oot ( oT 
S F I 
R ‘ R é 
0 2h I 20) 
i a) eo) ) 
Ry 2(Ro/Rp) Ry 2(R,/R)) 


45 


—— emma 
— 


= et 





The partial derivatives were evaluated numerically and were 


fmoeund to be: 


= Sones = = 1.375 
S F 
2h 

me onoeae oh 

oT EE ee 
F Ca 

2h 

= 0.01 
TO ED 


therefore: 


tu 
Het & or 


Aa 


il 


+0.51 


The mean velocity is given by: 


go water 1/2 
Gy = | 2 Ge ee 


aha 


The only parameter here that is open to uncertainty (due to 
the assumption that property values are correct) is the H, 
or height of water column in the micromanometer. Since 

QH = +0.02 inches: 


(oH 
Z H 


= +0.324 ft/s 
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ay 





since: 


VD 
ieee 
V 
= SM 
AN. = Neg CT) 
M 
= +635. 


Nusselt number is by definition: 


Nu = os 
therefore: 
ANu = Nu (42) 
=f (ant LAlyy 
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Deen LX OF 


Penne tOte—GORRECTOR EQUATION 


In order to simplify and shorten the process of adjusting 
strip voltages to obtain an isothermal surface as described 
in Chapter 3, a predictor-corrector program was developed and 
implemented on a TI-59 with PC-100C thermal printer. The pro- 
gram listing is included at the end of this Appendix. 

The program requires entry of the ambient temperature and 
desired temperatureof the strips. Then the strip number is 
entered and the program prompts for initial voltage drep and 
initial current of the strip. The program computes strip 


temperature using the formula: 


where C and S are taken from the calibration curves. 
Next the program calculates the voltage drop that would be 
expected at the desired strip temperature using: 


t/2 


E,I,(1,-T 


p? 
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Where Ry = required resistance, from calibration 


te 


5 


Fluid temperature 


Dlopemenmecalabration curve 


im”emvariac Om the test rig is then adjusted to produce 
this voltage drop across the strip. The strip must be allowed 
to come to thermal equilibrium after an adjustment is made. 

The calibration curves assume thermal equilibrium and are mean- 
ingless when temperature transients are taking place. 

The program listed requires 204 program storage locations 
and 76 memory storage registers. The TI-59 must be reparti- 
tioned to obtain adequate storage registers by entering 8 (2nd) 
meee, . this provides 319 program storage locations and 79 
data registers. All of the data listed must be entered to take 
Mmavantage of the full capability of the program. 

After the program and stored data are loaded the ambient 


temperature is specified by entering the temperature and pres- 


sing key "B". The desired strip temperature can then be 
entered and the key "C" depressed. The program is now ready 
for use. 


In use the number cf the desired strip is entered and key 
Mmeeuedepressed. The printer will then prompt "E?" and the initial 
Woltage drop should be entered and the "R/S" key depressed. The 
program will echo the current on the printer and will proceed 
to print the current temperature and desired voltage drop. The 


program will then pause and is ready for the next strip number. 
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FIGURE 2 - CROSS-SECTIONAL SKETCH OF THE MODEL 
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FIGURE 10 - TYPICAL VELOCITY WAVEFORMS 
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PIGURE 12 = SKETCH OF THE LUMPED RESISTANCE 
MODEL 
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